Directed Transport of Confined Brownian Particles with Torque 
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We investigate the influence of an additional torque on the motion of Brownian particles confined 
in a channel geometry with varying width. The particles are driven by random fluctuations modeled 
by an Ornstein-Uhlenbeck process (OUP) with given correlation time r c . The latter causes persistent 
motion and is implemented as (i) thermal noise in equilibrium and as (ii) noisy propulsion in 
nonequilibrium. In the nonthermal process a directed transport emerges, its properties are studied 
in detail with respect to the correlation time, the torque and the channel geometry. Eventually, the 
transport mechanism is traced back to a persistent sliding of particles along the even boundaries in 
contrast to scattered motion at uneven or rough ones. 

PACS numbers: 05.40.-a,87.17. Jj 



Introduction. Nowadays, ratchet models are em- 
ployed to explain a wide range of transportphenomena, 
in biophysics as well as in artificial devices 0, • Gener- 
ally, the term ratchet refers to nonequilibrium phenom- 
ena that can arise provided one of the space-time sym- 
metries that inhibits directed motion is broken 0. This 
can be caused either explicitly by an asymmetric, peri- 
odic structure or by an unbiased, but asymmetric drive 
Q. In our case, the symmetry breaking is realized by 
a nonvanishing mean torque together with the confining 
geometry. 

The torque leads to a circular motion and an alter- 
ation of the diffusive properties Q . Such motion appears 
in many real world phenomena. For example, it can be 
due to asymmetries in the propulsion of agents them- 
selves as occurring in the chemotaxis of sperm cells 
and nanorods 0,13. Also, a t orq ue appears for Janus 
particles under laser irradiation or due to an external 
magnetic field that is used to steer nanorods [l2j |. 

Furthermore, our particles are confined in an infinitely 
long channel with a periodically varying width, see Fig. 
[TJ One channel boundary is introduced as a reflecting 
disk. Later on it will be interchanged by a reflecting tri- 
angle, thereby adding another symmetry breaking. Fi- 
nally, we also consider a 'rough' lower wall where elastic 
scattering is modeled by equidistributcd reflection angles, 
regardless of the angle of incidence. 

Even though related studies of transport mechanisms 
possess some of the features, such as the channel layout 
in chaotic transport in Hamiltonian systems fl3l . Il4j , en- 
tropic particle transport [l5j or for confined biological 
agents [8}, its combination with a dissipative dynamics 
and a nonzero mean torque to realize directed transport 
is a novel feature. Hereby, the actual driving of the parti- 
cles is performed by the noise, which delivers a constant 
energy supply. 

Thermal Colored Noise. In our setting the par- 
ticle has unit mass m = 1, position r(t) = (x(t),y(t)) 
and velocity r(t) = v(t) = (v x (i),Vy(t)) at time t, that 
are driven by correlated noise e(t) = (e x (t), e y (t)). Such 
noise is called colored, its memory causes a persistent 
motion. We choose an OUP to drive our particles, which 
is a common implementation for it represents the most 




FIG. 1. Illustration of the channel geometry, a box of unit size 
whose left and right side are identified with each other, real- 
izing periodic boundary conditions. At the top and bottom 
specular reflections take place, (a) Lower boundary formed 
by a disk of radius R — 1.2 that protrudes by L = 0.2 into 
the box. (b) Edged boundary parameterized by H — 0.5 and 
W = 2/3. (c) Flat lower boundary with equidistributed scat- 
tering angles due to a 'rough' surface. 



natural continuous valued colored noise (l6|. In equi- 
librium a fluctuation dissipation relation (FDR) relates 
the autocorrelation function of the noise and the friction 
function 

te(a) ej -(a + *)) = ^e-l'l/^ = fcT 7 (|i|)% (1) 

with i,j £ {x,y}. Here the Einstein relation T>^ = jokT 
(70 - friction constant, k - Boltzmann constant, T - tem- 
perature) is used to relate the noise intensity T>£ to the 
heat bath from which the noise derives. Assuming that a 
FDR holds, we are forced to relinquish the Markov prop- 
erty by introducing a dissipative memory kernel 7(|t — s\) 
to govern the friction. We then speak of a generalized 
Langevin equation HH for the velocity, 



v(t) 



l(t-s)v(s)ds + e(t)+Clv(t). (2) 



The mean torque is realized by multiplication of v with 
the rotation matrix f2, resulting in the force flv = 
(Qv y ,—Qv x ). 

Without confinement, the stationary Fokker-Planck 
equation corresponding to this dynamics is easily solved. 
We then find that the velocity is Maxwell-Boltzmann dis- 
tributed. Numerically, we can confirm that this holds 
true in the channel geometry as well, both globally for 
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FIG. 2. (Color online) Trajectories of confined (cf. Fig. [T] (a)) particles exposed to nonthermal colored noise with several 
torques and correlation times at constant 70 = 0.2, T>^ — 0.04. The depicted time length is t\ = 33, the velocity is indicated by 
the arrows, between two consecutive arrows At = 1/4 has passed. 



the entire channel and in the vicinity of the boundaries. 
Hence we have a symmetrical velocity distribution and 
no directed transport occurs. 

The particles perform an undirected diffusive motion. 
As an estimate one can take the effective diffusion coef- 
ficient of thi s dy namics calculated for unbounded space, 
which reads [18| 



d cS = 2V ( 7o 2 + n 2 



(3) 



Apparently, it does not depend on the correlation time 
and coincides with the expression for particles driven by 
white noise. 

Nonthermal Colored Noise. In the second case we 
consider noisy propulsion in nonequilibrium, hence no 
FDR holds UM- A gain equipped with an OUP-driving, 
our system is governed by the set of Langcvin equations 



r(t) =v(t), v(t) = -j v{t) + n«(t) + e(t), 



e{t) = e(t) 




(4) 



Here denotes white Gaussian noise of zero mean 
with uncorrelated components, i.e. = and 

(Zi(sMs + t)) = 6(t)5 ij . 

In this way, we have effectively realized the motion 
of an active particle, it underlies a permanent corre- 
lated supply of energy from the OUP-noise which it dis- 
sipates during the persistent and curved motion due to 
Stokes friction. The persistence of this motion expresses 
nonequilibrium. With r c — > white noise follows realiz- 
ing again an equilibrium situation. 

By taking the Fourier transforms of cq. (TJ|, the spec- 
trum of the velocity can be traced back to that of the 
noise 



(u, 2 



7o 2 



n 2 ) 



-*S U6j (w). 



(5) 



( 7o 2 + n 2 - lo 2 ) 2 + 4w 2 7o 

With the Wiener-Khintchine theorem and the 
Lorentzian noise spectrum of the OUP S eiej (oj) = 

2T>£ (l + t 2 lo 2 ) 1 Sij, which follows from the left-hand 
side of eq. (Q]) , we obtain a rather lengthy expression for 
(v(s)v(s + t)). For t = it simplifies to 



(ViVj) 



7o (7ot c 



n 2 r 



1 u v 



(6) 



Obviously, the torque, the friction, and the correlation 
time dampen the particles' mean squared velocity. For a 
vanishing mean torque the result without external force 
field [lj| is reproduced. 

We notice that the unconfincd effective diffusion coef- 
ficient can be calculated from the Fourier transform of 
cq. (0) by integration over the time, which leads to the 
same result as for thermal colored noise (cf. eq. (|3J)). 

Directed Transport in Confined Geometries. 
Fig. [2] depicts some trajectories in the circular chan- 
nel for the nonthermal dynamics given by eq. (j4]) . Rising 
correlation times lead to a reduction of the total spatial 
displacements as do increasing mean torques, although 
less severely. With both influences acting on the parti- 
cles together, another effect can be observed: The parti- 
cles stay in the vicinity of the reflecting boundaries for a 
longer time and perform a curly hopping motion that 
changes to a narrow creeping for large r c [2fj|. Also, 
we notice that while the red trajectory (empty arrow- 
heads) is largely unbiased, with bigger torques the parti- 
cles travel a longer distance along the flat wall than along 
the curved wall. 

The velocity distributions for the same situation are 
depicted in Fig. [3] [2l|. While the velocity is Gaussian 
distributed for r c = 0, it remains symmetrical on an un- 
bounded plane or in the confined geometry without con- 
stant torque. In contrast, if all these influences act on 
the particles together, i.e. boundaries, finite correlation 
time and constant torque field, the velocity distribution 
loses its symmetry (green circles). For f2 = 1, the v x - 
distribution's left tail is lowered, while its right tail is 
raised. Thus, we have a net particle flux. Unlike for par- 
ticles in equilibrium, the velocity distributions are nar- 
rowing with growing |f2|. Furthermore, we notice that the 
confinement particularly leads to a narrowing in P(v y ), 
which now markedly differs from P{v x ). 

Let us now address the stationary particle fluxes J that 
pass in the x-direction through our channel. Herein, J is 
normalized with respect to the number of particles. 

The numerical results for the net transport are shown 
in Fig. 2] as a function of the correlation time for several 
values of the mean torque. The flux J has a maximum at 
approximately r c ps 2 — 3 for all depicted f2 values. For 
vanishing correlation times, J converges to zero, which 
is consistent in view of the white noise case that follows 
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FIG. 3. (Color online) Velocity distribution (v x - symbols, v y 
- lines) for particles driven by nonthermal colored noise with 
t c = 1, 70 = 0.2, T>£ = 0.04. Except for the red distribution 
(squares), particles are confined as in Fig. Q](a). For r c = 0, 
the same Gaussian distribution emerges regardless of fl. 
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FIG. 4. (Color online) Net transport for particles driven by 
nonthermal colored noise for several mean torques at 70 = 0.2, 
T>£ = 0.04. The channel geometry is specified in Fig. [l](a). 



in this limit. There, the equilibrium condition forbids 
directed transport as the noise strength in eq. (|4]) van- 
ishes. For large correlation times the net transport also 
converges to zero. 

The effect of the mean torque is examined in further 
detail for several geometries in Fig. [5] Evidently, J van- 
ishes for O = and also for strong torques regardless 
of the mean turning direction. The direction of the net 
transport changes with the orientation of f2. 

Given the emergence of a directed transport, we can 
justify its qualitative form, i.e. the asymptotics in r c and 
fi as well as the occurrence of a maximum, with regard to 
the behavior of the speed. As implied by eq. ^ it van- 
ishes both for large mean torques and correlation times. 
For transport to take place however, we need a torque to 
break the symmetry, and a nonzero correlation time for 
nonequilibrium. Hence we have no transport in either of 
the cases £7 — > or — > 00 and r c — > or r c — > 00, 
leaving only a maximum value in between, which has to 
change its direction with the sign if the sign of the torque 
switches. 

We now turn to the effects of the differing geometries in 
more depth. As we see in Fig. [5J the antisymmctrical be- 
havior in fl is conserved only as long as the lower bound- 
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FIG. 5. (Color online) Net transport for various confining 
geometries (cf. Fig. [TJ. All curves are at 70 = 0.2, T>$ = 0.04 
and t c = 3. 



ary keeps as left-right symmetry (cf. Fig. [T](a,c)). Oth- 
erwise, as for the edged boundaries (cf. Fig. Q] (b)), the 
heights of the positive and negative peaks differ, whereas 
the direction of the current coincides for W = 1/3 and 
W = 2/3 (i.e. is independent of the orientation of the 
sawteeth). Also, the position of the maximum moves to 
slightly smaller f2 values for the edged surface compared 
to the disk. 

If we narrow the channel by increasing the origin of the 
reflecting disk, the resulting transport curves also rise. 
However, their overall behavior remains; in particular the 
position of the maximum is barely altered and remains at 
Sl" lax ps 1 — 1.5. From this we conclude that the position 
of the fi max cannot be traced to some resonant-like effect 
between the channel width and the cyclotron radius R c = 
v/fl. Rather, the position of f2 max , as well as r" lax , is a 
tradeoff between the necessary symmetry breaking and 
the reduction of the particles' speed. 

Due to the unit scaling of the velocity, our net flux 
is identical to the mean velocity in the ^-direction, J = 
(v x ). Other correlation ratchets [22|, [23| display values 
for (v x ) in the same order of magnitude. 

Elucidation of the Transport Mechanism. There 
are three possible types of trajectories in our channel: 
particles in the middle perform a more or less circular 
motion without net bias, whose regularity is enhanced 
by ft and decreased by T>^ and 70. Particles near the 
top eventually collide with the wall and are reflected. A 
counter clockwise mean turning rate then results in a 
forward motion and multiple further reflections. If we 
consider a particle that is driven to the wall by a noise 
pointing into that particular direction, the noise works 
largely as a drag force after reflection. It is still oriented 
in the same direction, whereas the particle has changed 
its own. As a result, the velocity normal to the bound- 
ary decreases (cf. Fig. [3|). Once the noise orientation 
has changed (after about At ps t c ), the particle is driven 
away from the boundary, although this may be counter- 
acted by the particle's inertia together with the torque. 
Such situations are shown in Fig. [21 Clearly visible are 
the repeated reflections with ever decreasing radius and 
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speed at the top. 

Without additional obstacle, particles near the bottom 
perform an antagonistic motion that cancels out the for- 
ward flux. The obstacle on the other hand modulates the 
reflection angles depending on the point of incidence. As 
seen in Fig. [SJ the distance the particles have to travel 
until the next reflection varies depending on whether the 
surface tangential is increasing or decreasing. This con- 
stitutes a scattering effect that helps particles to escape 
from the vicinity of the boundary. 

Hence the backward flux is hindered and a positive net 
transport emerges. In this picture, we can also explain 
the augmentation of J with the narrowing of the channel 
(cf. Fig. [5]), for the particles cannot move as far without 
interaction with the boundaries. Thusly, the proportion 
of intermediate unbiased motion is reduced and the im- 
pact of the symmetry breaking is enhanced. 

We confirm this rationale by considering a box with 
a flat but rough lower boundary, modeled by random 
cquidistributed reflection angles (cf. Fig. Q] (c) and Fig. 
[5j black symbols). Obviously, this leads to a similar be- 
havior of J. 

Consequently, the flux is dominated by the sliding mo- 
tion along the flat boundary and affected by the per- 
sistence, the torque and the noise driven velocity. In 
contrast, the geometry of the confinement and the shape 
of the lower boundary appear to be less important and 
influence the total value of J rather than its form. To 
evaluate qualitatively the dynamical dependence of the 
current, we expect that J = J(r c , f2, |v|). It has to vanish 
if any of these three parameters vanishes, which favors a 
product of them. We can assume proportionality to the 



velocity, which we approximate by \J (vv) . Then, the de- 
pendence from the torque is met by an additional factor 
yJSi. The flux J rises linearly for small r c , in comparison 
to the numerical results this must be supplemented by 
an additional damping for large ones. Hence we propose 
J oc sgn(O) y|?2| (vv)t c cxp(— r c /4), which mimics the 
numerical results quite well. It yields a optimal torque 
at |f2 max | = t c 7 1 +7, where the transport has a maximum 
of J max cx t c exp(— r c /4). This behavior can also be seen 
in Fig. [U where f2 max wanders to slightly smaller values 
with r c , and in the height of the corresponding peaks. It 
can be interpreted as a time scale merging between the 
period due to the torque, which is proportional to I/O, 
and the two relaxation times t c and Ti ner tia = 1/7 that 
characterize the persistence of the motion. 

In conclusion, we have constructed a setup in which un- 
biased correlated noise is rectified in nonequilibrium and 
leads to the emergence of directed transport. To that 
end, we needed a nonzero mean torque and an asym- 
metric confining channel. We have given a qualitative 
explanation of the mechanism that leads to this sym- 
metry breaking, namely the particle hopping along the 
boundaries disrupted by the scattering effect of the re- 
flecting disk or triangle. Despite the general settings of 
the propulsion, particles and geometry, this investigation 
can help in the qualitative understanding and might be 
the starting point for more precise application oriented 
modeling of noisy active particles. 
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